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We study an instability of thin liquid- vapor layers bounded by rigid parallel walls 
from both below and above. In this system, the interfacial instability is induced by 
lateral vapor pressure fluctuation, which is in turn attributed to the effect of phase 
change: evaporation occurs at a hotter portion of the interface and condensation 
at a colder one. The high vapor pressure pushes the interface downward and the 
low one pulls it upward. A set of equations describing the temporal evolution of 
the interface of the liquid-vapor layers is derived. This model neglects the effect of 
mass loss or gain at the interface and guarantees the mass conservation of the liquid 
layer. The result of linear stability analysis of the model shows that the presence of 
the pressure dependence of the local saturation temperature mitigates the growth of 
long-wave disturbances. The thinner vapor layer enhances the vapor pressure effect. 
We find the stability criterion, which suggests that only slight temperature gradients 
are sufficient to overcome the gravitational effect for a water /vapor system. The 
same holds for the Rayleigh- Taylor unstable case, with a possibility that the vapor 
pressure effect may be weakened if the accommodation coefficient is below a certain 
critical value. 



1 



I. INTRODUCTION 



Thin liquid films have been intensively studied over the last decades. Many contributions 
have been devoted to them, owing to their technological importance and wide industrial 
applications. Their rich interfacial behaviors originate from combinations of various effects 
such as capillarity, intermolecular forces, thermocapillarity and gravity. One important effect 
among them is evaporation or condensation. It was often incorporated in the studies of thin 
films,-!^ and dewetting patterns resulting from drying of films were analyzed .-i^>^ 

Many past studies on evaporating or condensing liquid films have ignored the dynamics 
of the gas above the liquid film, assuming the infinitely deep gas phase.-i^i^ However, in this 
study we consider a situation where the gas phase is bounded by rigid parallel wall from 
above and has a finite depth comparable with the liquid one. The full linear stability analyses 
of this system were performed in several papers.-i^i^ii^iii Despite the apparent simplicity of 
the configuration, this system includes a free surface and interfacial boundary conditions 
involving phase change, and therefore is very sophisticated. In order to simplify this problem, 
we apply long-wave approximation to both layers. 

The advantage of the use of long-wave or lubrication approximation is the reduction of 
dimensionality: a one-dimensional (two-dimensional) film evolution equation can be derived 
in a two-dimensional (three-dimensional) system. Normally, only the dynamics of the liquid 
is considered, leading to a one-sided model. However, if the ambient gas layer is thin enough, 
a two-layer model would better describe the system. This was demonstrated by VanHook et 
al.,— who developed a two-layer theory to reproduce their experimental results. They showed 
that their two-layer model better predicts the onset of instability in their experiment than 
the corresponding one-layer model and also correctly describes the formation of localized 



2 



elevations. In their approach, only the heat conduction in the gas phase is taken into account, 
and the gas dynamics is ignored because the viscosity of the gas is much less than that of 
the liquid. Later, Merkt et al.— presented an evolution equation of the interface of two 
viscous fluid layers in the same geometry. Their model allows for the shear stress induced 
by the motion of the upper layer and therefore is reduced to the single layer equation in the 
limit of small viscosity of the upper layer. Although their goal is the observation of pattern 
formation in the long-time regime, the two-layer models have been applied to the cases of 
Rayleigh- Taylor instability^ii^ and ultrathin films.— 

Nevertheless, in the two-layer systems mentioned above there is no phase transformation 
at the interface. Here, we construct a two-layer theory for liquid-vapor layers which undergo 
phase change, using long-wave approximation. Note that the application of long-wave theory 
to the vapor phase was made in the study of film boiling.— If we take into account the effect 
of the mass flux across the interface, an instability peculiar to this system is expected, even 
for the presence of large disparity in viscosity and density between liquid and vapor; see 
Fig. [H The liquid film is initially in equilibrium with its vapor layer. If the liquid side is 
heated or the vapor side is cooled, evaporation occurs at a hotter portion of the interface and 
condensation at a colder one. Since the vapor layer is bounded, the vapor pressure becomes 
higher in the evaporating region and lower in the condensing one. According to this lateral 
vapor pressure gradient, the higher vapor pressure pushes the interface downward and the 
lower one pulls it upward. Then, the surface deflection is amplified. To our knowledge, 
this pressure-induced instability mechanism has not been considered in the past, because 
the uniform ambient vapor pressure has been assumed in previous studies of evaporating or 
condensing liquid films 

To derive the model, we require the interfacial boundary conditions such that the mass 
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FIG. 1: Instability mechanism of the present system. 

transfer occurs between the phases. We follow those of the earlier studies on evaporating 
or condensing liquid films^'^ except the thermodynamic relation at the interface. They used 
the linearized equation where the mass flux through the interface is proportional to the 
difference between the interfacial temperature and its saturation value corresponding to 
the surrounding vapor pressure, based on kinetic theory.— This relation cannot be directly 
applied to the present problem, because the local saturation temperature varies in the lateral 
direction, depending on the vapor pressure. Hence, we must modify the relation to take this 
effect into account. Fortunately, this can be easily done in the thermodynamic framework. 
For instance, Ajaev and Homsj»2S and Wayner— used a nonequilibrium thermodynamic 
relation including the saturation temperature variation due to capillarity and disjoining 
pressure. This effect was later included in the model of evaporating or condensing thin 
liquid films.— However, they assumed the vapor pressure to be constant. In our two-layer 
model, this relation should be extended in accordance with the vapor pressure variation. 
Thus, one of the purposes of this work is to investigate the effect of the vapor pressure 
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dependence of the local saturation temperature. Note that this effect was not considered 
in the film boiling case,— although the lateral vapor pressure variation, which drives the 
motion of the vapor, certainly exists in the boiling film. It is worthwhile noting that a 



23l . where the vapor concentration 



somewhat similar motivation to ours is found in Ref. 
in the ambient gas phase above the liquid film fluctuates and thereby the mass flux varies 
along the interface. However, they neglect the bulk gas dynamics itself and consider only 
the diffusion of the vapor. 

Here, we start with a more general nonequilibrium thermodynamic law, which reduces 
to, in the linear domain, a proportional connection between the interfacial mass flux and 
the difference of chemical potential in each phase.— From this law, we can naturally de- 
rive a thermodynamic relation similar to that of Ajaev and Homsy and Wayner. Moreover, 
the condition of local thermodynamic equilibrium, adopted in some papers on two-phase 
problem,->i2iiiii2^ is recovered by taking the appropriate limit of the derived relation. There- 
fore, the thermodynamic relation used here is also the extension of the interfacial equilibrium 
condition into nonequilibrium states. We note that more general formulation taking into 
account the nonequilibrium effect contains a temperature discontinuity at the liquid-vapor 
interface during evaporation or condensation, as was done in Ref. 2J. However, this tem- 
perature jump may be neglected unless the phase change occurs too rapidly.— 

In the derivation of our model, we manipulate the mass flux balance equation at the liquid- 
vapor interface. In the literature, the effect of mass loss or gain at the liquid surface due 
to evaporation or condensation has been included into the model through this equation by 
assuming that the vapor speed is much larger than the liquid one because the vapor density 
is much smaller than the liquid one.-i^ However, in this paper we show another interpretation 
of this equation as a consequence of order estimate. If we assume that the degree of the 



disparity in density between two phases is much greater than that in viscosity, which is vahd 
for most substances far from the critical point, we can find that the hquid velocity is not 
balanced by the vapor one in the mass balance equation. Instead, it is balanced by the 
interface velocity, which indicates that the effect of mass loss or gain can be neglected. This 
order estimate leads to the approximation of the mass balance equation, where the total 
mass of the liquid is conserved and the effect of the mass flux affects only the dynamics of 
the vapor. Under this approximation, we derive the model where the conservation of the 
total liquid mass is guaranteed, whereas the effect of evaporation or condensation remains 
in the vapor dynamics. 

The paper is organized as follows. The model is formulated in Sec. Ull where the interfacial 
boundary condition and scaling peculiar to this system mentioned above are introduced. 
Linear stability results are presented in Sec. Illlt including analyses of the Rayleigh- Taylor 
instability and the effect of degree of nonequilibrium on it. Section IIVI summarizes the 
results and future work. 

II. FORMULATION 

For simplicity we consider a two-dimensional system as in Fig. [1], where the horizontal 
bilayers, liquid and vapor of the same substance, are confined by rigid parallel walls from 
both below and above. We assume that the initial equilibrium thicknesses of the liquid and 
vapor layers, di and d^, are small enough to ignore the buoyancy effect. The temperatures 
of the liquid-side and vapor-side plates are controlled at Ti and T2. The x axis is taken to be 
parallel to the walls, and the z axis perpendicular to them. The z = plane corresponds to 
the boundary between the liquid and the liquid-side plate. The position of the liquid-vapor 
interface is described hj z = h{x,t). Gravity acts in the negative direction of the z axis. 
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A. Governing equations 

We assume that the equations of continuity for incompressible fluids and of the momen- 
tum and energy balance hold in each phase. They are given by, respectively, 

V ■ = 0, (la) 

p/3(5tv^ + V/3 ■ Vv^) = -Vpp + r]pV\f}- ppge^, (lb) 
dtT/3 + ■ VTp = KpV^Tp. (Ic) 

Here, = {up, wp), and Tp are velocity, pressure and temperature fields, respectively, in 
the (3 = {v, 1} phase, where v denotes the vapor and I the liquid. The differential operator 
is V = {dx,dz)- The coefficients p/j, 77^ and denote the density, dynamic viscosity and 
thermal diffusivity in the /3 phase, respectively, which are assumed to be constant in each 
phase. In Eq. ( llbp . g is the gravitational acceleration and the unit vector in the z 
direction. 

B. Boundary conditions 

At the walls (^ = and z = di + dg), we impose no-slip boundary conditions. Along with 
the temperature conditions prescribed above, they read 



V, = 0, Ti=Ti at z = 0, (2a) 

Vy = 0, Ty = T2 at z = di + dy. (2b) 

At the liquid- vapor interface z = h{x,t), the mass flux J must be conserved: 

J = Pvi^v ■ n - V/ • n) = pi{vi ■ n - V/ ■ n). (3) 
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Here, n is the unit normal vector directed toward the vapor, 

n = 4z£iiL, (4) 
and V/ represents the interface velocity, which satisfies the kinematic condition 

V/ ■ n = — , =■ (51 

We assume the continuity of the tangential velocity along the interface, 

■ t = V, ■ t, (6) 



where t is the unit tangent vector to the interface. 



The interfacial stress and energy balance equations read, respectively,-!^ 



(7) 



J(v, - V,) + {pi - p,)n - {2riiEi - 2t],E,) ■ n + 2aHn = 0, (8) 
J{L + |[(v, - V,) ■ n]2 - l[(v, - V,) ■ n]2} + A^VT^ ■ n - A.VT, ■ n 

+ [27]iEi ■ (vz - V,) - 2r],E, ■ (v, - v/)] • n = 0, (9) 

where E^, a, H, L and Xjs are the rate-of-strain tensor in the f3 phase, the surface tension, 
the mean curvature of the interface 

the latent heat and the thermal conductivity in the (3 phase, respectively. In Eq. ([8]), we 
ignore the thermocapillary (Marangoni) term and assume a as well as L and A/j to be con- 
stant for simplicity. The Marangoni effect on two-phase surfaces has often been neglected 
in the literature.-i^»^»i^ The recent investigations on the linearized systems of liquid-vapor 
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layers^iiSiiii and drop or bubble^^, show that the Marangoni effect has httle significance in 
pure two-phase coexisting states, because the hquid-vapor interface becomes almost isother- 
mal owing to the large entropy difference between the two phases. In the Appendix, we 
include the thermocapillary term in our model and examine its effect. We numerically find 
that the thermocapillarity makes little contribution to linear stability of a stationary state 
of the model at least in the physical situations considered here. The projection of Eq. ([8]) 
on the normal and tangent to the interface yields, respectively, 

-^(v; -Vy) ■n + pi-p.,-n - {2riiEi - 2r7^E„) ■ n + 2aH = 0, (11a) 
t ■ {2r]iEi - 2r],E,) ■ n = 0, (lib) 

where Eq. (Q was used in the second equation. Assuming the moderate phase change rate, 
the continuity of the temperature at the interface holds: 

Ti=T, = Tj. (12) 

Finally, in order to close the system we require an additional boundary condition, which 
relates to the interfacial thermodynamic state. In this study, we adopt the linearized phe- 
nomenological law such that the mass flux across the interface is proportional to deviation 
from local thermodynamic equilibrium:— 

J = K[ni{pi, Tj) - fiviPv, Tj)]. (13) 

Here, is the chemical potential in the (3 phase, which is a function of the pressure in 
the corresponding phase and the temperature at the interface. A proportionality coefficient 
K will be later specifled by analogy with the kinetic theory. We now expand the chemical 
potentials into Taylor series in this equation around their initial equilibrium value /xq with 
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respect to the variations of the pressure and the temperature, 

5pp = Pf3-Po, (14a) 
ST = Tj-T,at{Po), (14b) 

where Tsat{po) is the saturation temperature at the initial equihbrium pressure po- Using 
the Gibbs-Duhem relation for a one-component system, we obtain for each phase 

/i/3(p/3, Ti) = /io - S/sST + —6pf3, (15) 

where S/j is the entropy density of the /? phase. Substituting Eq. f|T5l) into Eq. f|T3l) yields 

J = K (as5T + -5pi - —5p^] , (16) 
V Pi Pv J 

where As = s^j — si is the entropy difference between the phases and related to the latent heat 
by L = Tsat{po)As. If we neglect the pressure terms on the right hand side of Eq. ( IT6|) . we 
recover the usual kinetic relation.->^'^*^>^ On the other hand, in the limit K ^ oo, Eq. ( |T6|l 
reduces to the condition of local thermodynamic equilibrium, used in several phase-boundary 
problems.-ii^iiii2^ Therefore, Eq. (fT6l) is an intermediate relation connecting the two different 
interfacial conditions appearing in the studies of two-phase systems with phase change. 

C. Dimensionless equations and parameters 



In order to nondimensionalize the above equations, we scale lengths, time, velocities, 

"IfL, JL, JL, AT and ^ 

r]i dipi dfpi diL 

where AT is the initial temperature difference across the liquid layer. We find AT together 



pressures, temperatures and mass flux by di, — — , - — , —J—, AT and — r^— , respectively. 



with TsatiPo), solving Eq. ffTcl) for both phases with the boundary conditions ([2]), ([9]) and 
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(fT2|) in the equilibrium steady state (J = 0), as follows: 

AT = Ti - r..,(p„) = j^(Ti - Tj). (17a) 
T.M = ^^fi^. (1Tb) 

Here, the dimensionless parameters A and d have been introduced. The definitions of di- 
mensionless parameters appearing in this paper are presented in Table [H Furthermore, we 
define the dimensionless pressure and temperature such that their initial equilibrium values 
at the interface, po and TgatiPo): correspond to in their new variables. In the following, we 
show the resulting nondimensionalized equations. 

First, the governing equations of the liquid layer ([1]) become 

V ■ = 0, (18a) 
dt^ri + V; ■ Vv; = -Vpi + V^Vi - Ge„ (18b) 
P{dtTi + vi-VTi) = V^Ti, (18c) 

and those of the vapor layer 

V-v^ = 0, (19a) 

p(5tv^ + • Vv^) = -Vp^ + r]V\^ - pGe,, (19b) 

P{dtT, + V, ■ VT,) = kV'T,. (19c) 

The boundary conditions at the walls ([2]) reduce to 

ui = wi = 0, Ti = l at z = 0, (20a) 

Uy = = 0, Ty = — Y at z = 1 + d, (20b) 

A 
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and those at the interface ([3D, (Q, ([11]), (Q, and ([IS]), respectively, 

EJ = p(v^ ■ n — V/ ■ n) = Vi ■ n — v/ ■ n, (21a) 

V, ■ n = ^ , (21b) 

■ t = • t, (21c) 

EJ{^^l - v„) ■ n + - - n ■ (2Ei - 2riE^) ■ n + 2SH = 0, (21d) 

t ■ (2E, - 2r]E^) ■ n = 0, (21e) 
J + ^{[K - V/) ■ n]2 - [(v, - V,) ■ n]^ + VT, ■ n - AVT, ■ n 

+ ^[E; • (vi - V,) - r]E, ■ (v, - vj)] ■ n = 0, (21f) 



Ti=T, = Tj, (21g) 



J = 



^ n / 1 

Ti + — [pi 



(21h) 



E V P 

In Eq. ( I21h[) we have introduced the dimensionless parameter i^, instead of K in Eq. ([16 



which has the dimension. The parameters K and K are related through K = diLK/Xi. 
The value of the parameter K defined in Table [T] is determined by the comparison with the 
kinetic theory (Hertz-Knudsen law). In its definition, a is the accommodation coefficient, 
m is the molecular mass of the fluid and ks is the Boltzmann constant. 

D. Long-wave asymptotics 

We apply the long- wave approximation to both layers.^ Letting a small parameter e be 
di/A, where A represents the characteristic lateral length scale, new space and time variables 
are introduced as 

x' = ex, z' = z, t' = et. (22) 
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TABLE I: Dimensionless parameters. 



Gravity 



G 



Density ratio 



Pi 



Liquid Prandtl number 



Evaporation number 



Surface tension 



E 



S 



n 



= ^ 
^ (ypidi 

rjf 



K = a 



{LdipiY 
PvdiL'^ 



m 



Dynamic viscosity ratio ij = — 

m 



Conductivity ratio A 



Diffusivity ratio 



Initial thickness ratio 



a7 



di 



hTsat V '^T^^bT. 



e 



sat 



AT 



sat 



This rescaUng indicates that the physical quantities vary much slower in the horizontal 
direction than in the vertical one. Assuming that e ^ 1, we expand the velocities, the 
pressures and the mass flux in powers of e as follows: 

ui = + (^Uh + • • • , M„ = e"^(M„o + eti^i + ■■■), 
wi = e{wi^ + ewi^ + ...), Wy = w^^ + ew^^ + . . . , 

(23) 

Pi = ^~\pio + mi + •••)> Pv = e"HPfo + ^1 + •••)> 

J = Jo + e Ji + . . . . 

Here, we required wp/up = (9(e) based on the continuity equations. To take the pressure 
effects into account, we chose pi,Pv = 0(e~^), where the pressures of both layers are taken 
as the same order in e because of the pressure balance. From this and the balance between 
the pressure and viscous dissipation terms in Eqs. fllSbl) and fll9bl) . |v;| and ?7|v^| must 
be the same order. Therefore, if we set i] = 0{e), |vi|/|v„| = 0{e) holds. In order for 
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all the physical effects to appear in the leading-order equations, some of the dimensionless 
parameters shown in Table [T] are scaled by e as 

G = e-^G, P = P, E = e^E, S = e'^S, U = e^Il, K = K, 

(24) 

p = e^p, r] = er], X = X, k = e~^K, d = d, 
where we set O ^ to neglect the molecular kinetic energy and viscous dissipation terms 
in Eq. (I21fj) . The tildes denote the quantities of order C(l), which will be used instead 
of the original dimensionless parameters. Here we assumed that the evaporation number 
E is of order O(e^), instead of 0{e) as in the work on evaporating or condensing liquid 
films.- This is because in our system only a small amount of evaporation or condensation 
is sufficient to drive the vapor dynamics owing to the very small vapor density compared 
with the liquid density (see a discussion below Eqs. (|27|) on decoupling of the mass balance 
equation fl21al) ). Since the density ratio p is much smaller than the dynamic viscosity one 
rj for most substances far from the critical point, we set p/rj = 0{e). Moreover, we took 
the thermal conductivity ratio A as 0{1) in the scalings fl2^ . although the real value of A 
is very small because in general the thermal conductivity of gas is much lower than that of 
liquid. We cannot take the limit A ^ because the temperature boundary condition at the 
vapor-side wall fl20bp contains the factor A^^ We substitute these scalings into the previous 



dimensionless equations and take the limit e — 0, so that only the leading-order terms in e 
are left in the equations. Hereafter, we shall omit the primes, the tildes and the subscripts 
0, unless otherwise stated. 
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The leading-order governing equations are 



dxUi + d^wi = 0, (25a) 

d.pi = dluu (25b) 

d,vi + G = 0, (25c) 

dlTi = 0, (25d) 



for the hquid layer {0 < z < h) from Eqs. (ITSl) . and 

dxUv + dzWv = 0, (26a) 

dxPv = vdzUy, (26b) 

d,p, = 0, (26c) 

d^T, = 0, (26d) 

for the vapor layer {h < z < 1 + d) from Eqs. (flQl) . Whereas the wall boundary conditions 
(ppj) remain unchanged, those at the interface (z = h) from Eqs. (12 ip result in 

dth = -uidxh + wi, (27a) 

EJ = p{-u^dxh + Wy), (27b) 

= 0, (27c) 

p,-p, + ^92/i = 0, (27d) 

d^ui = rjd^u^ (27e) 

J + d^Ti - \d,T, = 0, (27f ) 



Ti=T, = Tj, (27g) 

n 



J = K(Ti- 4tP, ) . (27h) 
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Note that from Eq. (12 lap with Eq. (121bp we obtain the two decoupled equations fl27ap and 



(]27bl) . because the hquid and vapor velocities in Eq. (12 lap do not balance under the scalings 
(^^. In other words, since we have scaled the density ratio as O(e^), the terms representing 
mass loss or gain become the next order in e in Eq. fl27ap . and hence are discarded there. 
The next order mass balance gives Eq. fl27bp . where the effect of the interface velocity in 



Eq. fl21al) has been neglected compared to that of the vapor velocity. Therefore, Eq. fl27ap 
implies that the total mass of the liquid layer is conserved at the leading order in e, while 
the effect of the mass flux affects only the dynamics of the vapor layer through Eq. fl27bp . 
In addition, the vapor recoil term in the normal stress balance (121 dp and the liquid pressure 
term in the thermodynamic relation (121hp have disappeared with these scalings. 



The origin of the decoupling of the mass flux balance equation (]21ap is more specifically 
explained as follows. As is mentioned below Eq. (I23p . |v;|/|v^| is the same order as the 
dynamic viscosity ratio rj. From the fact that p <^ rj, it follows that p ^ |v/|/|vy| or 
pWv\ ^ |v;|. For the second equality of Eq. (1 2 lap to be true, the hquid velocity v/ ■ n 
must be balanced by the interface velocity v/ ■ n, which is represented by Eq. (I27ap through 
Eq. (I21bp . In Eq. fl27bp . v/ ■ n has disappeared because from Eq. (I27ap it is the same order 



as ■ n, much smaller than v„ ■ n from the small viscosity ratio 77 -C 1. Therefore, in order 
for the decoupling of Eq. (12 lap into Eqs. (I27ap and (]27bp to be valid it is essential that 
p <^ 7] <^ 1 he the case. 



Solving Eqs. (]25dp and (]26dp with the boundary conditions (I2U]) and ( |27gP yields the 



temperature gradients in both layers 



a.r,^Zi^. ..T„^-IL1^. (28) 
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Substituting these equations into Eq. (127f[) and eliminating J using Eq. (127hp . we obtain 



, ^ n \ Ti-l XTj + d 



Then, the surface temperature T/ can be exphcitly expressed as 



(29) 



Ti 



1 



1 + Kh 



\h 



'l + d){h-l) KYI 
l + d-h ^ pE^\ 



(30) 



l + d-h 

Substituting this equation into Eq. (]27hp again, we find the expression for the mass fiux 



J 



K 



l + Kh + 



Xh 



[l + d){h-l) 
l + d-h 



1 



\h 



l + d-h J pE 



n 



(31) 



l + d-h 

From Eqs. fl25cp and fl26cp . we can find that the horizontal pressure gradients dxPi and 
dxPv do not depend on the vertical coordinate. Then, we can twice integrate Eqs. ( ]25bl) and 
fl26bp in the z direction. Using the boundary conditions fl2Ul) . fl27cl) and fl27ep . we obtain 



ui = ^dxPiz^ + ciz, r]Uy = ^dxPv{l + d - zf + C2(l + d - z) 



(32) 



with 



ci{x,t) = -^{1 + d- h)dxPv - hdxPi, C2{x,t) = ^{l + d - h)dxPv (33) 

The expressions for the vertical velocities wi and rjw^ immediately follow from the integration 
of the continuity equations fl25ap and fl26ap with the no-slip boundary conditions fl20l) : 

wi = -Idlpiz' - l-dxCiz^ riw, = ldlp,il + d-zY-l-d.M^ + d-zY. (34) 

Z / 



Integration of Eqs. fl2^ and fM) with Eq. (1275]) gives the relation between the liquid and 



vapor pressure 



dxPi = dxPv - Sdlh + Gdxh 



(35) 
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Substituting Eqs. §^ and with Eq. ([33]) into Eqs. ([27a]) and fl27b]) finally yields a set 
of equations, respectively, 



dih = 



+ -^da:{Gh - Sd^h) 



(36) 



EJ = -^d,[{l + d-hfd^p,l (37) 



where the liquid pressure gradient dxPi has been eliminated using Eq. fl35l) . Equations (|3T]) 
and (l37|l can be combined to eliminate J: 

\h 

Equations (I5B1) and (1351) compose a closed system for the unknown variables h{x,t) and 
Pt;(x, t). The first term in square brackets of Eq. ([36]) describes the effect of the lateral vapor 
pressure gradient, while the second that of gravity and the surface tension. The second 
term on the left hand side of Eq. fl38|) represents that of the variation of the local saturation 
temperature due to the vapor pressure fluctuation. Equation ([5^ is written in the conserved 
form for h, because we have neglected the effect of mass loss or gain by decoupling the mass 
flux balance equation (12 lap as before. 

III. LINEAR STABILITY ANALYSIS 



The set of Eqs. f [36]) and ([38]) has a stationary solution h = 1 and = 0. We perturb 
this state by 

h{x,t) = 1 + h exp{ikx + ut) , (39a) 

Pv{x,t) = p exp{ikx + ut) , (39b) 
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where h and p are infinitesimal quantities. Linearizing the system gives the following growth 
rate: 

Ah'^ 1 

. = ^-iP(G + ^P), (40) 

with 

a 

Here, positive (negative) values of uj indicate instability growth (decay). The parameters A 
and fcg represent the effect of lateral vapor pressure fluctuation and that of local saturation 
temperature variation, respectively. From the dispersion relation ( HOl) . one can easily find 
that the growth rate vanishes in the limit — > 0. Nevertheless, if fcg were not present 
in Eq. fHOj) . the finite growth rate would remain for k ^ 0. Therefore, the presence of 
the saturation temperature variation mitigates long-wave growth rates. Specifically, the 
saturation temperature is increased (decreased) in higher (lower) vapor pressure regions and 
thereby the rate of evaporation or condensation is reduced. This effect is prominent for 
long-wave disturbances. Notice that our model does not admit a quasisteady solution of 
fiat moving interface even if ko = 0, because the lateral uniformity leads to dth = from 
Eq. fl36|) : this is a direct consequence of the neglect of mass loss or gain. We can consider 
only the long- wave limit k ^ 0, where a; 7^ if /cq = 0. 

A. Superheated or supercooled state 

To quantify the above results, we consider the water/vapor system at 100 °C and 1 atm. 
Using the material properties shown in Table [Tll we plot the growth rates (140|) in Fig. [2l 
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where we set di = 10~^ m, AT = 0.1 °C, g = 9.8 m/s^ and a = 1. In the experiment, it is 
the temperature difference between the plates Ti — T2, not across the hquid layer AT that 
can be controlled. However, in the following we fix AT so that the vertical temperature 
gradient in each phase is constant when we vary the value of the initial thickness ratio d. 
From Eq. (I17ap AT = 0.1 °C corresponds to T1—T2 = 2.8 °C for = 1 in this system. In the 
short-wave regime {k ^ 1), the growth rate is negative because of the effect of the surface 
tension {S in Eq. (HOj) ). whereas it is reduced in the long- wave regime {k <^ 1) owing to that 
of the saturation temperature dependence on the vapor pressure [k^ or 11 in Eqs. fHUl) and 
(]41bl) ). as mentioned above. 

TABLE II: Physical properties of water/vapor at 100 °C and 1 atm, identical with Table I of Ref.|8. 

pi = 960 kg/m^ = 0.6 kg/m^ L = 2.3 x 10^ J/kg 

r]i = 2.9 X 10"^ kg/m s rj^ = 1.3 x lO'^ kg/m s cr = 5.8 x 10"^ N/m 

A, = 6.8 X 10-1 J/m s °C A^, = 2.5 x lO'^ J/m s °C 

Ki = 1.7 X 10-'^ mVs K„ = 2.0 X 10"^ mVs 



The three dispersion curves with different values of d in Fig. [2] suggest that the instability 
is enhanced for the thinner vapor layer, which reflects the fact that A in the dispersion 
relation ( HOl) is a monotonically decreasing function of d according to Eq. (14 lap . As can be 
seen from Eq. ( I40l) . the factor A roughly represents the intensity of the instability unless k^ 
is very large. Physically, the friction at the wall of the vapor side prevents the vapor flow 
from mitigating the lateral vapor pressure gradient, and this effect is more pronounced for 
the narrower vapor layer. Therefore, the destabilizing effect of the lateral vapor pressure 
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FIG. 2: Growth rates a; = ea) versus wavenumber k = ek for the water/vapor system at 100 °C 
and 1 atm, where uo and k correspond to uo and k in the dispersion relation (|40|) . Here di = 10~^m, 
AT = 0.1 °C, g = 9.8 m/s^, a = 1 and d = 0.5, 1 and 2 from top to bottom. 

variation is stronger for the thinner vapor layer. The stabihzing role of the vapor flow 
is also understandable by considering the viscosity of the vapor. From Eqs. (HI]), A and 
also /cq are proportional to the viscosity ratio r], which implies that increasing the vapor 
viscosity intensifies the instability. The large viscosity of the vapor weakens the vapor flow 
from the equation of the vapor motion ( ]26bl) and also from the tangential stress boundary 
condition fl27ep . Thus, the reduction of the destabilizing vapor pressure effect by the vapor 
flow is suppressed. Although both decreasing d and increasing rj impede the vapor flow, 
the dependence of A or /cq on them is different, because the vapor flow is prevented by the 
different mechanism. 

The above argument is confirmed by directly estimating the vapor flow. Integration of the 
horizontal component of the vapor velocity (1521) in the vertical direction gives the expression 
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for the total vapor flux 



qv= u^dz = -T7p(l + d- hfd^py. (42) 

Jh ^^'n 

Similarly, the hquid flux is obtained as 

qi = r Uidz = _ /^ + 3(l + ^) ^2g^^^ _ }^f^3g^(^Q^ _ gQ2j^y (43) 

In terms of and g;, the model equations (136!) and (l37Il are rewritten as 

EJ 

dth = -d^qi, = d^q^. (44) 

P 

If we regard h as nearly 1, the vapor and liquid fluxes are 

qv ^ "^"^^P^' * - ~^2^ d^Pv - ^d^iGh - Sdlh). (45) 

Substituting Eqs. (H5l) into Eqs. ( 14^ . one can find the origin of the intensity of the instability 
A: the factor r]/d^ in the expression for A fl41al) arises from the vapor flux and A + 3d from 
the liquid flux. The rest is a contribution from phase change or the temperature difference. 
Thus, it is shown that the stabilizing vapor flow is responsible for the strong dependence of 
the instability intensity on the initial thickness ratio and also for that on the viscosity ratio. 

It is worthwhile noting that the role of the vapor flow presented here is different from 
that described in the study of full linear stability analysis of a similar bilayer system, where 
the penetration of the fluid at the walls is allowed.^ The authors of Ref. |8| attributed the 
stabilizing effect of the vapor flow to the convection of heat: the vapor flow convects heat 
from a hotter portion of the interface to a colder one. Their description cannot be applied to 
our system, because the term representing heat convection does not appear in the leading- 
order equations, (I25dp and fl26dl) . of the long- wave approximation. 
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From Eq. (HOl) . we can determine the cutoff wavenumber kc analytically. By setting u = 0, 
we obtain 



1 + 2- 



'iki/2 + by 



(46) 



where 



«=2^' ^=2^- 

If the right hand side of Eq. fB6l) does not take a real positive value, the cutoff wavenumber 
kc does not exist and hence the dispersion curve never crosses the line u = 0. This indicates 
that the system is linearly stable, because from Eq. ( HOj) the surface tension makes the 
growth rate negative in the short-wave limit k ^ oo. Since k^ is negative when a — bk^ < 0, 
the criterion for the linear stability is expressed as 

a ?>A {A + ?,d){l + d) pE ' ^ ^ 

This inequality is equivalent to the condition that the growth rate is negative at infinitesimal 
wavenumber, which can be found by Taylor expansion of Eq. (HUj) around A; = 0. This 
stability condition is independent of the factors representing the vapor flux contribution to 
A, rj/d^, and the effect of degree of nonequilibrium, ^ _^ ^ _^ because /cq also has the 
same factors and they are canceled out. For the same condition as before, the left hand side 
of Eq. ( I48l) becomes 7.3 x 10~^ for d = 1. From this evaluation, it is concluded that only 
slight temperature gradients are sufficient to overcome the stabilizing gravitational effect 
in the realistic system. Nevertheless, it is noted that the unstable modes of infinitesimal 
wavenumbers can be eliminated by appropriately modulating the horizontal scale of the 
system. Therefore, the stability criterion for such a system should be weaker than Eq. (jUj). 
The expression for the fastest growing mode kmax can be also analytically obtained from 
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Eq. (HOj) . Straightforward calculation yields 

1/3 



2kl + h \{kl-hf 1 2 1 akl[A{kl-hf + 27akl 



^max — o + \ 27 2*^^° ~'~ 6 



1/3 



+ ^(M^ + i„t„^-l./^M!l(M^^)!±^^ . (49) 

27 2 6V 3 

Unlike the results of linear stability analysis of the other film equations, a simple relation 
between kc and kmax cannot be established for our model. Equation (H9|) seems to be too 
complicated to find any asymptotic form. However, in the special case /cq = ^) reduces to 
a simple form 

kl,, = -b + [aby/'. (50) 

This special case is realizable for the water/vapor system considered if we set di = 2.3 x 10~^ 
m for d = 1 and K — > oo. 

B. Rayleigh- Taylor instability 

We also investigate the Rayleigh- Taylor instability of the system,-i2iii2i'i^i^ii^ii^ where grav- 
ity acts toward the vapor side. The system of interest is the case that the layers are heated 
from below or cooled from above, so that the stabilizing effect of evaporation or condensa- 
tion counteracts the destabilizing one of gravity. In this subsection, we consider the balance 
between these two effects by changing the signs of gravity and the temperature difference. 

Before starting the analysis, we show the difference from the model of evaporating or 
condensing liquid films with infinitely deep vapor layer.- In Ref. |6|, the authors made the 
two assumptions: the much larger gas layer depth than the liquid one and the neglect of the 
latent heat in the temperature boundary condition at the interface ([9]) (Eqs. (3) or (4) in 
Ref. |6|). To compare their dispersion relation with ours, we abandon these assumptions and 
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show their dispersion relation derived without imposing them. In our notation, it reads 



\ABM\ + ^k^iG~Se), 



(51) 



with 



\A 



1 + d 



K 



BM\ 



d 



\E\ 



(52) 



The assumptions which they made are equivalent to ^ 1 and K 1. Hence, by taking 
these limit we obtain = K\E\, which is consistent with the result of Ref. [g. However, 

for the water/vapor system considered here K = 1.2 x lO^a for di = 10~^m. Therefore their 
assumption of the negligible latent heat {K ^ 1) is questionable unless the accommodation 
coefficient a is very small. Their dispersion relation (!5T!) has the same form as ours (l40l) . 
except for the presence of k^. However, their definition of A in Eq. (132]) is different from 
ours (141 aj) because of the different mechanism: local mass loss or gain at the interface is 
the main stability mechanism in Ref. [g. In contrast, this effect has been neglected in our 

system compared to that of the vapor pressure fluctuation. This is also conflrmed by taking 

\Abm\ P d^ 



the ratio between both A: 



Note that this ratio is independent of K. 



\A\ T]4: + 3d 

For the water /vapor system, I^ba/ |/|A| ~ 0.011 for d = 2. Therefore |Abm| ^ |^| holds in 
the system considered here. Here the condition for the neglect of the effect of mass loss or 
gain, p <^ 1], is required, as was mentioned before. However, if d is much larger, the effect 
of mass loss or gain can be comparable to that of vapor pressure fluctuation, because the 
latter effect is much more weakened for a thicker vapor layer. 
The cutoff wavenumber for this case is given by 



^2 



\b\ 



-1 



1 + 2 



— |a| + \b\k\ 



{ki/2-\b\y 



1 ± Wl + 2 



-\a\ + 

m - kl/2Y 



for |>|6|, 



for |<|6|. 



(53) 
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The second line of this equation suggests the existence of two cutoff wavenumbers when 
— \a\ + |6|A;g < 0. Recall that —\a\ + |6|A;q < is the condition that the growth rate is 
negative at infinitesimal wavenumber. Therefore, the instability starts around k = for the 
first case of Eq. fl53p . whereas finite modes between the two cutoff wavenumbers are unstable 
for the second. The finite critical wavenumber where the instability starts for the second 
case is kcrit = \/\b\ - kl/2. 

In seeking the critical condition for the stability, it is desirable to vary the liquid depth 
di independently. In the dispersion relation fHU]) . we have three dimensionless parameters 
depending on di, k"^, G and S. Since S is proportional to di, we choose S" as a control 
parameter and the remaining two parameters are scaled by S to obtain new parameters 
independent of di: 

^-^3' ii-^ii, fco-^^O- ^2^4 l + K+X/d^^- ^^^^ 

For simplicity, here we ignore the effect of K by assuming the local thermodynamic equi- 

K 

librium ^ oo or — = 1 in Eqs. fHTl) because K ^ 1 for a = 1 as mentioned 

1 + K + X/d ^ 

above. 

The stability condition for the Rayleigh- Taylor instability is obtained from the cutoff 
wavenumber (153!) similarly to the previous case and expressed in terms of the above new 
parameters as 



4(rf + A) 


G*|n*5 


(4 + 3rf)(l + rf) 


p 


E 





< 1 for ^ < ^, 



CI 



(55) 



I, r]{A + ?,d){l + d) 12\E\S^ 
where = {k"^* /\G*\y^\ Here the two cases of this condition each correspond to those 

of the cutoff wavenumber fl53|) . The first line represents the condition that fc^ is negative, 

which is essentially identical with the previous stability criterion (HHjl . while the second the 
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one that fc^ is not real. In Fig. [3] we plot the neutral stability curves in S vs. |E| plane 
with the remaining parameters fixed as shown in Table IIIII One can see the defiections of 
the curves, which correspond to the transitional points {S = Sc) between the two criteria 
in Eq. (1551) . The stability curve for d = 2 passes near the point of di = 1.0 x lO^'^m and 
|AT| = 0.1 °C. For this point, the system is stable for both d = 0.5 and d = 1. From 
Fig. [3] the stable region is wider for the thinner vapor layer, suggesting the enhancement of 
the stabilizing effect of lateral vapor pressure fluctuation. Figure H] displays the dispersion 
curves for the Ray leigh- Taylor unstable case with d = 2 and |AT| = 0.1 °C, \g\ = 9.8 m/s^, 
a = 1 around di = 1.0 x 10~^m. The fastest growing mode can be obtained from Eq. (149|) 
with a,b < 0. 

10-^ 

X 10-2 

tq 

~ 10-^ 

10-^ 

10"^ 

0.1 1 
S X 10"^ 

FIG. 3: Stability diagram in 5 = e~"^S vs. \E\ = e^lE] plane for the Rayleigh- Taylor unstable case. 
Here d = 0.5, 1 and 2 from bottom to top. The vertical line corresponds to di = 1.0 x 10~^m and 
the horizontal one I ATI = 0.1 °C. 



iATi = o.rc 




stable 


di=1.0xlO"V 

unstable 



27 



TABLE III: Values of the dimensionless parameters for the water/vapor system at 100 °C and 1 
atm under terrestrial conditions. Here, the tildes are restored to avoid confusion. 




1 2 3 4 5 

kx 10^ 

FIG. 4: Growth rates to versus wavenumber k for the Rayleigh- Taylor unstable case. Here d = 2, 
\AT\ = 0.1 °C, 1^1 = 9.8 m/s^ a = 1 and = 0.9 x lO^^m, 1.0 x 10-^m and 1.1 x lO^^m from 
bottom to top. 

C. Effect of degree of nonequilibrium on the Rayleigh- Taylor instability 

In Fig. [3], we did not consider the effect of nonequilibrium because K ^ 1 for a = 1 in 
the system considered. However, the accommodation coefficient a can be much less than 
unity and thereby K might approach one. Here, we examine its effect on the stability 
of the system. As was mentioned above, the degree of nonequilibrium K does not enter 
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the stability condition (HHj) for the thermodynamic unstable case. For the Rayleigh- Taylor 
unstable case, however, the criterion (155!) includes K only in the lower line {S > Sc) when 
K is finite. Then, it becomes 





1 + K+ X/d{\G*\S^ + k^*f 


77(4 + 3d)(l + rf) 


K 12\E\S^ 



< 1 for ^ > Sr. 



(56) 



Therefore, the coefficient of the left hand side of this inequality increases as K decreases, 
which may render the system unstable even if it is stable at large K. In Fig. [5] we show 
the stability diagram in a vs. |E| plane for the water /vapor system. For this system, the 
boundaries of the stability are almost constant if a > 10^^. The steep changes of the neutral 
stability curves are found for a < 10^^, where the system becomes unstable for the same 
value of the evaporation number E. The physical meaning of this behavior is that as the 
resistance to evaporation or condensation is increased the stabilizing vapor pressure effect 
no longer overcomes the destabilizing gravitational effect. 

2 



1.5 



o 



X 1 

tq 



0.5 



Stable 




unstable 



10- 



10" 



10"^ 
a 



10" 



10" 



FIG. 5: Dependence of the stability on the accommodation coefficient. Here d = 2 and di = 
0.8 X 10~^m, 0.9 X 10~^m and 1.0 x 10^'^m from bottom to top. The horizontal line indicates 
lATl = 0.1 °C. 
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IV. CONCLUSION 



We have discussed the instabihty of thin hquid-vapor layers bounded by rigid parallel 
walls from both below and above. In this system, the interfacial instability is induced by 
lateral vapor pressure fluctuation, which is in turn attributed to the effect of phase change: 
the vapor pressure becomes higher at an evaporating portion of the interface and vice versa. 
The liquid is driven away from the higher pressure place and pulled up to the lower pressure 
one. This pressure-induced instability mechanism has not been considered in the past. 

In the formulation, the interfacial boundary condition taking into account the pressure 
dependence of the local saturation temperature, Eq. f|T6|) . was imposed. The relation f|T6]) 
is an extension of the kinetic law and the condition of local thermodynamic equilibrium, 
conventionally used in the literature. We applied the long-wave approximation to both 
liquid and vapor layers, assuming that the layers are too thin for thermal convection to 
occur. The choice of scalings we adopted (l24l) allows us to decouple the mass flux balance 
equation ([3]), resulting in the neglect of mass loss or gain through the interface by evaporation 
or condensation at the leading order. In the dimensional form, we can formally write the 
decoupled equations fl27al) and (]27bl) as 



where the first and second equations correspond to the leading-order and next-order equa- 
tions in e. The condition for this approximation isp^r^^l: p rj for the neglect of 
the effect of mass loss or gain and 77 <^ 1 for that of the interface velocity in the second 
equation of Eq. (IFTI) . To our knowledge, this decoupling approximation of the mass flux 
balance between two phases has never been considered and might have a possibility to be 




(57) 
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applied to other phase-change problems where p <^ri holds. 

As a result, a set of the equations describing the temporal evolution of the interface of the 
liquid-vapor layers has been derived. One of the equations fl36p is written in the conserved 
form for the film thickness, as is expected from the decoupling of the mass flux balance 
mentioned above. Therefore, the total mass of the liquid layer is conserved in our model. 
On the other hand, the effect of the lateral vapor pressure gradient induced by phase change 
is included in Eq. ( I36i) and the vapor pressure is enslaved to the film thickness through 
Eq. fl55]) . The result of the linear stability analysis of this model shows that the presence of 
local saturation temperature variation by the vapor pressure mitigates the growth of long- 
wave disturbances. The instability is enhanced for the smaller initial thickness ratio d and 
larger dynamic viscosity ratio r], which increase the kinetic resistance to the vapor flow. The 
role of the vapor flow is to mitigate the effect of the lateral vapor pressure variation, which 

n 

is different from that described in Ref. 8. We also determined the criterion for the linear 
stability of the system and found that only slight temperature gradients are sufficient to 
overcome the stabilizing gravitational effect for the water/vapor system. 

We also considered the Ray leigh- Taylor instability of the system. Here, the stabilizing 
vapor pressure effect is balanced with the destabilizing gravitational effect under very small 
temperature difference between the plates. Again, the thinner vapor layer strengthens the 
stabilizing effect of lateral vapor pressure fluctuation. However, for this case the instability 
domain may be widened if the accommodation coefficient a is below a certain critical value. 
This value is about 10~^ for the water/vapor system of di = 10~^m. 

In this paper, we addressed only the linear stability of our model. The next step is to 
proceed to the nonlinear analysis of the model. We shall investigate the behavior of the 
solution of the equations in the nonlinear regime by means of numerical simulation. Three 
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dimensional computation of the model will reveal the possibility of occurrence of pattern 
formation as reported in Ref . |6 . 

The validity of the long-wave approximation, on which our model is based, should be 
examined. To this aim, the comparison with the full numerical simulation is required. 
In particular, we do not know the critical condition for the onset of small-scale cellular 
convection in liquid-vapor layers, which is neglected within the framework of the long-wave 
approximation. There should exist a critical temperature gradient or thicknesses of the 
layers (critical Rayleigh number) for the transition between the conductive and convective 
states of the temperature fields if the buoyancy effect is taken into consideration. It would 
be also of interest to make a comparison with the existing full linear stability analyses^i^ and 
even weakly nonlinear analysis^"^ of bilayer systems. The effect of heat convection should be 
estimated to clarify the role of the vapor flow, which is different between theirs and ours. 

Finally, we make two remarks on our model. First, we set the scalings in e on the 
dimensionless material parameters as Eq. fl2^ assuming their values for the water/vapor 
system at 100 °C and 1 atm as a representative substance. Hence, if the material properties 
are considerably changed (e.g. near the critical point), we must reset the scalings appropriate 
for the relevant values, which may lead to different evolution equations. Second, we can 
incorporate the effect of intermolecular forces as disjoining pressure in the formulation. 
This effect will be dominant for layers of thicknesses below lOOnm, as in Ref. y and llTl . 
However, for this scale the application of continuum theory to the vapor layer would not 
be valid, because the mean free path of a gas molecule amounts to about 60 or 70nm at 
atmospheric pressure and becomes much larger at reduced pressure. 
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APPENDIX: EFFECT OF THERMOCAPILLARITY 

In the main text, we ignore the thermocapillarity for simplicity. Here we examine its 
effect on our model. If the thermocapillarity is present, we must add the thermocapillary 
term in the stress balance at the interface ([8]): 

J{y^i - V.) + (pi - p,)n - {27]iEi - 2r/,E,) ■ n + 2aHn + (V,a)t = 0, (58) 

where V<; = t ■ V. Accordingly, Eq. (jllbp is modified as 

t ■ {2r]iEi - 2r/,E,) ■ n = V,a. (59) 

Here we use the linear approximation for the dependence of the surface tension on the 
temperature, 

o- = cTq - 7(T7 - TsatiPo)), (60) 

where ctq is the surface tension at the initial equihbrium temperature Tsat{po), corresponding 

to a in the main text. The coefficient 7 = — ; — is positive for most common substances. 

dTj 
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Then, the nondimensionahzed tangential stress balance equation fl21ep becomes 



t ■ - 27]E,) ■ n = -MV.Tj, 



(61) 



where the Marangoni number is defined by 



M 



(62) 



In order to include the thermocapillary effect in the model, the Marangoni number is scaled 



as 



M = e-^M, 



under the long-wave approximation. Then, Eq. fl27ep reduces to 



(63) 



(64) 



where the tilde over M is omitted. We note here that the temperature dependence of 
the surface tension (!60!) also allows the tangential variation of the normal capillary stress 
in Eq. ( |58l) . However, this effect is safely neglected in the framework of the long- wave 
approximation (see Ref. [jj). Using Eq. one of the coefficients in Eq. (155]) is modified as 

1 



ci{x,t) 



[1 + d- h)d^p, - hd.,pi - Md^Ti, 



(65) 



while the other does not change. As a result, only one of the equations 0361) contains the 
Marangoni term: 



/z + 3(l + ci) ^,^^^^ ^ ^^^^^^ _ ^ ^Md^Tj 



dth = 

Here, the interfacial temperature gradient is calculated from Eq. (1301) as 

1 



(66) 



1 + Kh + 
dxh 



Xh 



l + d)dd,h KU^ . 

H wOx{hp^) 



1 + d- h 



1 + Kh + 



Xh 



1 + d-h 



K + 



1 + d-h)^ 
1 + d 



[1 + d-hy 



pE 
:X 



(67) 



:i+d){h-i) Ku 

H —hp^ 



1 + d-h 



pE 
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Note that if we set = 0, Eq. ([67]) is 



{1 + d-h + Xhy ^ ^ 



The Marangoni term with this temperature gradient in Eq. fl66|) is identical with that of the 
two-layer model without phase change in Ref. Il2. 



Under Eqs. ( l66|l and (138|1 . the dispersion relation (l40l) is altered as 

d 



where 



^,^^{l + d) K 



6MU K 

(4 + U)E + 



P 1 + K + 



A 



(2 L d J 



(70) 



Here only the terms in the first line of Eq. ( l67l) contribute to the Marangoni terms in Eqs. ( !69|l 
and fITU]) . The Marangoni term in Eq. flB^ arises from the surface deflection, whereas that 
in Eq. ( ITOi) from the saturation temperature variation due to the vapor pressure gradient. 
If the Marangoni number M vanishes, Eqs. (1691) and (!70!) are identical to Eqs. (140!) and 
(I41al) . We have Eq. (!69|) be the same form as Eq. (HOj) by introducing a new dimensionless 
parameter 

G' = G-5i±i?^L_. (71) 

We estimate the relative importance of the Marangoni effect on A' and G' by numerically 
comparing the two terms in Eqs. (!70|) and (!7T!) . For the water/vapor system at 100 °C and 1 

atm, 7 = 2 X 10"^ N/m°C, so that M = 23e for di = IQ-^m and AT = 0.1 °C. Therefore, it 

6Mn K .„,„4 31 + dM 1 

follows that — — = 4.7 x 10"^ and -- 



p{A + 3d)El + K + X/d 1 + K + X 2 d Gl + K + X/d 

0.62 for d = 1. If ^ 1, the thermocapillary effect seems to be negligible in 

1 + K + X 



the linear regime, which agrees with the former results on the two-phase problems.- 



Yet both values increase as di decreases or AT increases because the former is proportional 
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to (independent of AT) and the latter to d'['^ and AT. If K decreases, the former 
decreases and the latter increases. However, the latter does not diverge as does the left 
hand side of Eq. fl56l) . In both thermodynamic and Rayleigh- Taylor unstable cases, where 
the vapor pressure and gravity effects counteract each other, the Marangoni effect acts as 
amplifying the former and diminishing the latter. Finally, we note that we do not know 
whether the Marangoni effect is negligible in the nonlinear regime, which will be ascertained 
in the ongoing numerical analysis. 
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